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We discuss the chiral phase transition in hot and dense QCD with three light flavors. Inspired 
by the well-known fact that the £/a(1) anomaly could induce first order phase transitions, we study 
the effect of the possible restoration of the Ua(1) symmetry at finite density. In particular, we 
explore the link between the Ua(1) restoration and the recent lattice QCD results of de Forcrand 
and Philipsen, in which the first order phase transition region near zero chemical potential (fi) 
shrinks in the quark mass and fi space when n is increased. Starting from the Ginzburg-Landau 
theory for general discussions, we then use the Nambu-Jona-Lasinio model for quantitative studies. 
With the partial Ua(1) restoration modeled by the density dependent 't Hooft interaction, we fit 
the shrinking of the first order region found in de Forcrand and Philipsen's lattice calculation at 
low fi. At higher /x, the first order region might shrink or expand, depending on the scenarios. 
This raises the possibility that despite the shrinking of the first order region at lower fi, the QCD 
critical end point might still exist due to the expansion at higher fi. In this case, very high precision 
lattice data will be needed to detect the recently observed back-bending of the critical surface with 
the currently available analytic continuation or Taylor expansion approaches. Lattice computations 
could, however, test whether the (7a(1) restoration is responsible for the shrinking of the first order 
region by computing the r\ mass or the topological susceptibility at small [i. 

PACS numbers: 12.38.Aw,11.10.Wx,11.30.Rd,12.38.Gc 



I. INTRODUCTION 



The J7a(1) anomaly is an interesting phenomenon of 
Quantum Chromodynamics (QCD). In the chiral limit 
with three massless quark flavors, QCD has the chiral 
SUl(3) ® SUji(3) symmetry. This symmetry is sponta- 
neously broken in the QCD vacuum, giving rise to eight 
massless Nambu-Goldstone bosons. The J7a(1) symme- 
try, on the other hand, is broken not spontaneously but 
explicitly due to the quantum anomaly. Thus, the rf me- 
son is not a Nambu-Goldstone boson. It remains massive 
in the chiral limit. This is the so-called £7a(1) problem 
[l| and its consequence 't Hooft showed that instan- 
tons, which are topological configurations of the classical 
gluon field, are related to the VpSX) anomaly. He also 
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constructed an effective quark interaction which breaks 
the E/a(1) symmetry Q (see also [H). The coupling con- 
stant of this (Kobayashi-Maskawa) 't Hooft interaction 
measures the effective strength of the UpSX) anomaly rel- 
evant to the hadron spectrum. 

It is interesting to investigate the effective restoration 
of the U pSX) symmetry at finite temperature (T) and/or 
quark chemical potential (/i) induced by the decrease of 
instantons 0, Q , even though the triangle anomaly is in- 
dependent of the infrared scale T @, [1] • Consequences 
of the E/a(1) restoration have been investigated if the 
restoration is associated with the chiral transition [{J [l?} • 
In Ref. Q, based on the instanton liquid model, it was 
speculated that drastic L/a(1) restoration at the chiral 
transition could be expected because the chiral transition 
is caused by the rearrangement of the instanton configu- 
rations in that model. 

An analysis using the Nambu-Jona-Lasinio (NJL) 
model was given in Ref. |iq |. where the topological sus- 
ceptibility xt, a correlator of topological charges, was 
considered. In the large N c (number of colors) limit, xt 
is related to the rj' mass through the Witten-Veneziano 
mass formula [11], 2NfXt/f% = m l + m l> ~ 2™^, so 
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it can be used to probe the J7a(1) anomaly. The NJL 
model calculation |lOj | reproduced the lattice data [IH 
above the critical temperature up to 1.5 times the chiral 
phase transition temperature with temperature indepen- 
dent 't Hooft coupling constant. This implies that, at 
least in the NJL model, the effective E/a(1) restoration 
does not necessarily take place near the chiral transition. 

At finite density, the medium response of the £/a(1) 
anomaly is even less understood, largely because of the 
lack of lattice data due to the fermion sign problem 
[T3 |. Reliable analyses are available only at asymptot- 
ically high density where QCD is perturbative. In Refs. 
[3 111, the i mass in the Color-Flavor Locked (CFL) 
phase was calculated and shown to decrease as the den- 
sity increases. This corresponds to the realization of ef- 
fective Ua(1) restoration. There are efforts to test this 
restoration hypothesis by pro bing the decrease in the rj 
and rj' masses 0, GJ, Gj, [H [HT H3] in medium. It is 
claimed that the decrease can be studied experimentally 
through observing the pos sible formation of -q- and rj- 
mesic nuclei [H M, HI ll, E3, HI . 

It is interesting to speculate how the anomalous vio- 
lation and restoration of J7a(1) symmetry could change 
the QCD phase diagram (see Ref. (29[ for example). In 
this work, we will not discuss the possible changes of 
phases and phase boundaries in the intermediate den- 
sity. Instead, we will focus on the issue regarding the 
existence of the QCD critical end point (CEP). In the 
conventional three flavor picture [30, HH HH , there is a 
first order boundary of the chiral phase transition sep- 
arating the hadronic and quark phases starting from a 
point with zero T but non-zero fi to a point (/i c ,T c ), 
then the QCD phase transition becomes a crossover at 
lower \x. This end point of the first order phase transi- 
tion, (/i c , T c ), is called the CEP. This CEP has an second 
order phase transition despite the finite quark masses. 
The search of the QCD CEP is a priority in the next 
phase of the RHIC running. 

However, this conventional picture of the QCD phase 
diagram and the existence of the CEP are challenged 
by the recent lattice QCD results of de Forcrand and 
Philipsen [33j . The assertion in Ref. [33j can be explained 
by the so-called Columbia plot, where each point in the 
parameter space of quark masses and fi is marked by 
its order of phase transition when T is increased. We 
will study QCD with three light flavors, u, d and s, with 
isospin symmetry m u = md = TO u d- When fi = 0, the 
chiral limit point (to u( j = Tn s = 0) has a first order phase 
transition when T is increased due to symmetry reasons 
Q. Around the chiral limit, there exists a finite area 
such that each point in this area represents a first order 
phase transition. The boundary of this area is marked by 
a critical curve. Each point on this critical curve is of a 
second order phase transition. Beyond the critical curve, 
the theory has a crossover until all the quarks become 
heavy such that the theory is close to a pure Yang-Mills 
theory and phase transitions can then take place again. 
In the following, we will just focus on the light quark 



region. 

At /I = 0, the physical point of quark masses is lo- 
cated in the crossover region, meaning that QCD has a 
crossover at fi = 0. If the CEP exists at finite fi, then the 
physical point should enter the first order phase transi- 
tion region at finite \x. However, in Ref. (33|, it was found 
that near \i = 0, the region of the first order phase tran- 
sition shrinks as fi increases. Note that this computation 
is not directly carried out with finite fi QCD due to the 
fermion sign problem. Instead, methods with analyti- 
cal continuation from imaginary fi and with derivatives 
computed at fi = are used. Thus, essentially it is the 
curvature at fi = that was computed. Although that 
result of Ref. (HI disfavors the existence of the CEP, to 
give a definite answer to whether the CEP exists or not, 
computations at higher fi are necessary (note that other 
lattice results using different approaches are consistent 
with the existence of the CEP ll, [H, [H, [13 ) • 

Given that lattice computations at higher fi are still 
challenging, it is our hope that model calculations might 
shed light on this problem. 

Inspired by the fact that the J7a(1) anomaly could in- 
duce the first order transition in massless three-flavor 
QCD, we study the effect of the possible restoration of the 
Ua(1) symmetry at finite density. Starting from general 
discussions using the Ginzburg-Landau theory, we then 
use the NJL model for quantitative studies. With the 
partial J7a(1) restoration modeled by the density depen- 
dent 't Hooft interaction, we fit the negative curvature of 
the critical surface at fi = obtained in Ref. [H[ . Finally, 
we discuss the behavior at higher fi within this model. 



II. GINZBURG-LANDAU THEORY ANALYSIS 

In this section, we use the Ginzburg-Landau (GL) the- 
ory to demonstrate the special role of the E/a(1) anomaly 
in chiral phase transitions, following the framework of 
Refs. d^HIillil- 

The effective potential ttcL for the order parameter 
(chiral condensate) field $y = (<7j(l — 75)^) in QCD 
with three flavors (Nf = 3) is parametrized by a set of 
operators satisfying the SU h (3)®SU R (3)®Uv(l)®Z(N { ) 
symmetry, where Z(N{) is the remaining symmetry of 
J7a(1) after it is broken by quantum anomaly. Here, we 
will use the mean field approximation and neglect the 
space time dependence of $. This approximation neglects 
soft-mode fluctuations, which are large near second or- 
der phase transitions, but it is still useful to explore the 
phase structure. Also, we will Taylor expand Qgl in an- 
alytic functions of the order parameter. The following 
analysis helps us to understand our NJL result in later 
discussions. We have 

Q GL = ^Tr$t$+^(Tr$t$) 2 + ^Tr($t$) 2 

~ (det* + det$ f ) - ^Trh ($ + $ f ) , (1) 



3 



where the last term breaks chiral symmetry explicitly 
with ho oc diag(m u , mj, m s ). Adding higher dimensional 
operators does not change the analysis qualitatively. The 
determinant term simulates the Ua(1) anomaly and has 
the Z(N{) symmetry. We have Co > at T = to yield 
a finite rf mass in the chiral limit. 

For simplicity we restrict ourselves to the flavor SU (3) 
symmetric case: m u = = m s (that is, on the diagonal 
in the Columbia plot). Then we have <I> = diag(<r, ct, ct). 
The GL functional is reduced to 



3 / 9 3 \ 

^gl = 2 ao<j2 ~ Ca(j3 + [ 41 bl + 41 &2 ) u ~~ 3/loCr 

= ^ct 2 - ^cct 3 + J&ct 4 - /jct , (2) 

where b > such that the free energy is bounded from 
below. 

(i) Chiral limit (h = 0) 

We first consider the chiral limit case with h — 0, where 
interesting results can be obtained. 

If a > and c 7^ at the phase transition tempera- 
ture, then Eq. ^ has two local minima and the phase 
transition is of first order. The first order transition per- 
sists against an external field h until it is washed out 
at sufficiently large h. Thus, the first order region has 
some finite extent from rn u( j = m s = in the Columbia 
plot. Going back to h — 0, the two minima are located 
at cti = and 02 = (c + y/c 2 — Aab)/2b. The chiral con- 
densate is 



ACT = CT2 — CTl 



Aab 



2b 



(3) 



Note that as c — > 0, if there is a phase transition, then 
a — > at the phase transition point and the phase tran- 
sition becomes second order. Here we have assumed as 
c — ► 0, b is still positive and that the higher order terms in 
ct can further be neglected. Otherwise first order phase 
transition is still possible through the inclusion of the 
higher order terms such as ct 6 . Only in this simple case, 
the absence of the J7a(1) anomaly leads to the disappear- 
ance of the first order transition. 

For simplicity let us assume an (unphysical) extreme 
case that the fj, dependence lies only in the anomaly 
term c and c — > at /1 = (J, les . Then as we go from 
fi = to higher fi on the fj, axis, the strength of the 
first order transition gets weakened. At fi = /i ros , the 
first order transition disappears completely, and above 
fi rC s, the transition turns into second order. However, if 
fluctuations are taken into account, the above mean-field 
picture is modified. The second order transition above 
ju res has the SU L (3) <g> SU R {3) ® U A (l) symmetry. The 
renormalization group analysis tells us that there is no 
infrared stable fixed point for this universality class Q. 
This means that the fluctuations wash out the critical 
point and make it a fluctuation induced first order tran- 
sition. Thus, at m u d = m s = 0, the transition remains 



first order on the /i axis even if the £7a(1) symmetry is 
restored. This implies that as we go to ^ r cs, the first or- 
der region around m u d = m s = would ever shrink, but 
the critical surface can never touch the fj, axis but only 
approaches it asymptotically. 

The above discussion gives a clear picture about the 
eventual shrinking of the first order region if ?Za(1) is 
completely restored above a certain fi (if terms higher 
order than Eq. Q can be neglected) . For partial restora- 
tion with the GL coefficients depending on [i, the even- 
tual shrinking does not have to happen. This can be seen 
from Act in Eq. ([3]). If c decreases, Act can still increase 
if the parameter b (or a) decreases sufficiently fast. This 
means that even if the Ua(1) symmetry is partially re- 
stored when /i increases, the first order region may still 
expand. In fact, this is confirmed in the NJL calculation 
in the next section. In contrast, even without the E/a(1) 
symmetry restoration, the critical surface can still shrink 
if b or a increases due to some mechanism. This is similar 
to what happens if a repulsive interaction between vector 
currents is added to the NJL model, as demonstrated in 
Ref. [4l[. Thus, the Ua(1) restoration is just one possible 
explanation to the shrinkage of the first order region at 
small /i. 

(ii) Finite current quark mass (h 7^ 0) 

The different scenarios mentioned above can also be 
seen in the following analysis involving finite quark 
masses. At the critical end point, the GL functional takes 
the form of 

Q GL = -b (ct - CT ) 4 + d 

= -b (ct 4 - 4ct ct 3 + 6ct 2 ct 2 - 4ct^ct + ct 4 ) + d(4) 

Comparing this with Eq. ([2]), we have 

13 1 1 

b(T ° = 3 C ' 2 ba ° = 2°' ba ° = h ' 4 &CT ° + d = - ( 5 ) 

These are the conditions to determine the critical sur- 
face. However, there are six parameters — a, 6, c, h, <tq 
and d, while there are only four relations among them 
in Eq. ([5]). Thus, as mentioned above, one cannot iden- 
tify what causes the shrinking of the first order region 
without extra inputs. 



III. NJL MODEL ANALYSIS 

In this section, we give a quantitative analysis using the 
NJL model [H, IH, [ij] with the partial La(1) restoration 
modeled by the density dependent 't Hooft interaction. 
We will fit the curvature of the shrinkage of the first order 
region at fi — to the lattice QCD result of Ref. [33| . 
We then discuss the behaviors at higher fi. 
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A. Model Setting 



The NJL Lagrangian is given by 



C 

Co 



Cq + £4 + 

q{vy-d-rh)q, 



a 

3d [det%(l - 75)^ +h.c] . 



(6) 
(7) 

(8) 

(9) 



The kinetic term £0 includes the current quark mass ma- 
trix rh which breaks chiral symmetry explicitly. The sum 
of the two four Fermion contact interaction terms in £4 
is chirally symmetric. It gives the attractive interaction 
responsible for the spontaneous chiral symmetry break- 
ing when its strength exceeds a certain critical value. In 
this minimal model, the other types of four Fermion con- 
tact interactions of the same mass dimension are not in- 
cluded because they are heavier excitations below the 
phase transition. The six Fermion determinant term is 
the 't Hooft interaction which breaks the J7a(1) symme- 
try. 

The theory has five parameters: m uc j, m s , g$, <7 D , and 
the cut-off A. They fit the following physical quantities in 
the vacuum: m n = 138 MeV, /„. = 93 MeV, m K = 495.7 
MeV, and m v , = 957.5 MeV with ra ud = 5.5 MeV fixed. 
The resulting values are m s = 135.7 MeV, gs = 3.67A~ 2 , 
and A = 631.4 MeV [Z^. Note that in this paper, the 
't Hooft coupling constant is determined by fitting the 
T]' mass, g D = — 9.29A~ 5 = g&, where we have defined 
5d as the value of g D at the vacuum. Alternatively, one 
can choose to fit it to the mass boundary of the critical 
surface using lattice results. 

We will follow the standard procedure to calculate the 
effective potential of the NJL model in the mean-field 
approximation [3(| 52, 51, 53 • From the effective poten- 
tial one can construct the critical surface. The technical 
details will be omitted here. 

To apply the NJL model to finite T and /x, we will 
include the fi dependence in <? D , g D = <7 D (/i), while the 
other parameters m u( j, rn s , gs, and A remains indepen- 
dent of T and [i. The T dependence of c/ D is small in 
this model and is thus neglected. Note that in Ref. [nj, 
the T independent g D was found to be able to repro- 
duce the T dependence of the topological susceptibility. 
Without knowing the functional form of g D (/i), we plot 
the critical curves in the (//, m = m u d = m s ) space with 
different constant g D 's in Fig. [1] (the Columbia plot with 
3d (A*) = 3d is shown in Fig. 3 of Ref. [Hj].) Each curve 
ends with a point beyond which first order phase tran- 
sition ceases to exist, but crossover transition can still 
happen. This plot shows that the first order regions ex- 
pand as /i increases and larger the value of g D , larger the 
first order region. One immediately realizes, if <7 D is a 
monotonically decreasing function of fx, then the critical 
curve will bend to the left (corresponding to the shrink- 
ing of the first order region) at small ix, as the lattice 




0.5 1 1.5 2 2.5 3 3.5 4 

Quark mass m ud -m s [MeV] 



FIG. 1: The critical curves of constant with g^, — 
gu, 0.8(?d, 0.4<?d, and 0.2go, respectively. If is a mono- 
tonically decreasing function of /i, then the critical curve will 
bend to the left (corresponding to shrinking of the first order 
region) at small g,, as the lattice result of [33| has shown. At 
higher g, whether the curve can keep bending to the left or 
eventually bend to the right (corresponding to shrinking or 
eventual expansion of the first order region), depends on the 
functional form of <?d(a*)- 



result of [33| has shown. But, at higher /i, whether the 
critical surface could keep shrinking or rather expands, 
depends on the functional form of g D (/i), i.e., how fast 
the J7a(1) symmetry gets restored. 

As a special example, we will use the ansatz, 



<7D(M) = <?D(0)e-^° , 



(10) 



with one free parameter fj,o to parametrize the £/a(1) 
symmetry restoration. This form is motivated by the 
Gaussian suppression of the instanton density due to De- 
bye screening [H 53, H 51, HO] ■ 



B. Lattice Data 

Before going to the numerical results of the NJL model, 
we discuss the lattice data of the critical surface [331 ] . 
The critical curve, which is the intersect of the SU(3) 
symmetric plane, m — m n< \ = m s , to the critical surface 
can be Taylor expanded near li = 0, in even powers due 
to CP symmetry: 



m c (fi) 
m c (0) 



= 1-3.3(3) 



JL 

ttT c 



-47(20) 



JL 

ttT c 



(11) 



It is not clear what the radius of convergence of the above 
expansion is. By requiring the third term to be smaller 
than the second term (which are both ~ 23% of the first 
term), we obtain fi < 90-100 MeV. 
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The negative signs for the fj, 2 and /x 4 terms suggest the 
shrinking of the first order region at small /x. The sign of 
the /x 6 term was also asserted to be negative, making the 
shrinking even more serious, although the value of the 
coefficient has not yet been determined [33] . However, as 
we discussed above, we do not expect the critical surface 
to touch the /x axis at zero quark masses. Thus, there 
must be terms at higher order in the /x expansion with 
positive prefactors which change the critical surface from 
the charecteristic front bending to back bending at some 
density. 



C. Numerical Results in the NJL Model 

Now let us discuss the numerical results of the NJL 
model. In Fig. [21 we show the result of the critical curve 
in the (/x, m = m U A_= m s ) plane for various values of the 
free parameter /xo [51( . The shaded region represents the 

/j [MeV] 
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0.5 


1 1.5 2 2.5 



m c (/x) /m c (0) 



FIG. 2: The critical curves in the (/x,m = m u( j = m s ) plane 
with <?d (0) = <?d for several values of /xo. The shaded region 
represents the lattice result of Eq. (|lip . The lattice value of 
T c depends on /x. At /i = 0, T c — 135 MeV which gives a 
typical value for T c . 

lattice result of Eq. (fTTj) . The band shows the errors in 
the expansion coefficients but not from the higher order 
terms or the systematic errors. 

Several curves are plotted in the figure. As we go from 
the left-hand curve to the right-hand one, /xo increases 
from 31.6 MeV to infinity. The infinite /iq corresponds 
to the density independent <?q. The critical curves termi- 
nate at the critical chemical potentials at T = 0, above 
which first order phase transition ceases to exist, but 
crossover transition can still happen. We see that the 
values of /xq can be divided into three regions according 



to the qualitative behaviors of the corresponding critical 
curves. 

(i) The small (1q region: In this region, the suppression 
of the ?Za(1) anomaly is so strong that the first order 
region shrinks quite rapidly. The surface keeps shrinking 
and approaches the /x axis asymptotically. 

(ii) The intermediate /xo region: At low density, the 
first order region shrinks and then expands again at high 
density. This back-bending behavior is due to the den- 
sity effect. In general, stronger first order transition is 
favored at higher density, as discussed in detail in Ref. 
PH . In the language of the GL theory, this means that 
the higher order terms enter the thermodynamic poten- 
tial. Thus at finite density, there is a competition be- 
tween the J7a(1) restoration which makes c smaller and 
the density effect which generates higher order terms. In 
this intermediate /io region, the density effect overcomes 
the Ux(l) restoration at high density, resulting in the 
back-bending structure. 

(iii) The large /x region: When /x is large, the U&(1) 
restoration is weak. Thus, the first order region does not 
shrink and instead expands monotonically. This corre- 
sponds to the conventional scenario with constant in 
which the CEP could exist. 

Now let us compare the curves with the lattice data. 
We see that the curves with (j,q = 253 ~ 379 MeV de- 
scribes the lattice data well for /x < 70 MeV. It is no- 
table that this value of /xo is close to a rough estimate 
by no ~ l/(\/A 7 fPo) ~ 380 MeV, where N t = 3 and 
Po ~ 0.3 fm is the typical instanton size. For 70 MeV 
< fj, < 100 MeV, /x = 189 - 253 MeV describes the 
lattice data better. 

In our calculation m c (0) ~ 1.1 MeV is obtained while 
the lattice data gives m c (0) ~ 14 MeV (extracted from 
Fig. 9 of the first paper in Ref. [33j). One can argue 
that since m n ' is heavier than the cutoff A of the theory, 
it should not be used to fit the parameters of the the- 
ory. Instead, one can use the m c (0) computed in lattice 
QCD. To explore the effects of this new set of parame- 
ters, we use (/{j (0) = 3.5(/d and keep the other param- 
eters the same. This parameter set gives m c (0) ~ 26 
MeV, not quite the same as the lattice to c (0), but rather 
the Columbia plot at /x = is very similar to that of 
[33l ] . In Fig. [31 we plot the corresponding critical curves. 
The qualitative features are the same as those in Fig. [2] 
There is a range of no (253 ~ 379 MeV) where the curva- 
ture from lattice QCD is reproduced. In the range below 
/io = 253 MeV, the first order region does not expand 
substantially while in the range above (1q = 379 MeV, 
it expands beyond the critical mass at /1 = 0. Also, the 
results of Figs. [H and O suggest that the range of /i 
that corresponds to the curvature from lattice QCD is 
insensitive to the size of <?q (0) or m c (0) in our model. 

To further explore the existence of the CEP within 
this model, one needs to go beyond the m u d = m s limit. 
Again, using <?q (0) = 3.5<?d, we show the critical surface 
in the (/x, m u d, m a ) space with /ig constrained by lattice 
data. In the panel (a) of Fig. H £t = 253 MeV, the 
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FIG. 3: Same as Fig. [2] but with (0) = 3.5g D . 



C/a(1) restoration is strong enough such that the criti- 
cal surface does not intersect the physical quark mass 
line which is denoted by the thick dash line in the plot. 
Thus, CEP does not exist in this case. In the panel (b) 
(/xq = 379 MeV), on the other hand, the critical surface 
intersects the physical quark mass line. Thus, the CEP 
exists. 

A few comments are in order: 

(a) The most interesting feature of our result is the 
back-bending behavior of the critical surface. In our 
model we use [Iq to parameterize the E/a(1) restoration 
and fit it to lattice data. But similar behavior was also 
seen in another model with the repulsive vector-vector 
current interaction added to the NJL model [4l|. Thus, 
the back-bending seems to be a generic feature of the 
NJL model, and it could even be a property of QCD. 
This raises the possibility that despite the shrinking of 
the first order region at lower /i, the QCD critical end 
point might still exist due to the expansion at higher /i. 
In this case, it might be challenging for lattice QCD cal- 
culations using the analytic continuation and the Taylor 
expansion to probe the back-bending behavior. It will 
be worthwhile to investigate how precise the lattice com- 
putation should be in order to detect or rule out the 
back-bending scenario. 

(b) As mentioned above, the reason for shrinking of 
the first order region is uncertain yet. Both partial £/a(1) 
restoration and the vector-vector repulsive four-fermion 
interaction can give this effect. Lattice computations 
could, however, test whether the ?Za(1) restoration is 
responsible for the shrinking by computing m,y or the 
topological susceptibility at small /i. 




FIG. 4: The critical surface in the (ju, m„d, m s ) space with 
the 't Hooft coupling g^ (0) = 3.5<?d to better describe the 
critical line in lattice at /i = 0. The panel (a) is for fio = 253 
MeV and the panel (b) is for ^ = 379 MeV. The physical 
quark mass lines are denoted by the thick dashed lines. The 
physical quark mass line intersects the critical surface which 
means the CEP exists in the panel (b). In contrast, CEP 
does not exist in the panel (a). In the panel (a), the critical 
surface has a boundary curve above which we have no phase 
transition. It looks like a parallel curve but actually it has a 
small but finite dependence on fi. 



IV. SUMMARY 

We have discussed the chiral phase transition in hot 
and dense QCD with three light flavors. Inspired by the 
well-known fact that the JZa(1) anomaly could induce 
first order phase transitions, we have studied the effect of 
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the possible restoration of the c7a(1) symmetry at finite 
density. In particular, we explored the link between the 
C/a(1) restoration and the recent lattice QCD results of 
de Forcrand and Philipsen, in which the first order phase 
transition region near zero chemical potential (fi) shrinks 
in the quark mass and /i space when fj, is increased. Start- 
ing from the Ginzburg-Landau theory for general discus- 
sions, we then used the Nambu-Jona-Lasinio model for 
quantitative studies. With the partial J7a(1) restoration 
modeled by the density dependent 't Hooft interaction, 
we have fit the shrinking of the first order region found 
in de Forcrand and Philipsen's lattice calculation at low 
fi. At higher fi, the first order region might shrink or 
expand, depending on the scenarios. This raises the pos- 
sibility that despite the shrinking of the first order re- 
gion at lower //, the QCD critical end point might still 
exist due to the expansion at higher /i. In this case, very 
high precision lattice data will be needed to detect the 
back-bending of the critical surface with the currently 
available analytic continuation or Taylor expansion ap- 



proaches. Finally, since the rf mass and the topological 
susceptibility are sensitive to the strength of the /7a (1) 
anomaly, lattice computations of these quantities at small 
fi could check whether the strength of the E/a(1) anomaly 
is reduced when \i is increased. These calculations, how- 
ever, have to be carried out using either imaginary chem- 
ical potential or derivative expansions due to the fermion 
sign problem with finite \x. 
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